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Abstract 

During  the  past  years  we  have  intensively  studied  basic  properties  and  field  solutions 
in  homogeneous  bianisotropic  media.  We  started  from  simple  isotropic  media  to  ever  more 
general  bianisotropic  media.  This  study  has  led  to  some  classification  of  bianisotropic  media, 
which  resembles  a pedigree.  The  pedigree  contains  two  branches.  The  first  branch  are  so- 
called  self-dual  media,  these  are  generalisations  of  the  chiral  medium.  The  second  branch 
axe  factorizable  media,  these  are  generalisations  of  the  uniaxial  anisotropic  media.  We  think 
to  have  reached  some  consensus  with  respect  to  the  pedigree,  i.e.  we  think  to  have  found  the 
most  general  media  in  each  of  the  two  branches. 

1.  Introduction 

In  this  contribution  we  want  to  report  on  some  of  the  findings  we  obtained  during  the  past  eight 
years.  During  this  period  we  have  been  investigating  homogeneous  bianisotropic  media  with 
constitutive  relations  of  the  form 

D = e E + £ H,  B = C-E  + /r-H,  (1) 

where  e,  /i,  £ and  £ are  the  medium  dyadics. 

We  were  interested  in  finding  basic  field  solutions  in  these  media.  The  fields  for  any  given 
source  can  be  found  by  integration  from  the  fields  of  an  elementary  dipole  source,  i.e.  from  the 
Green  dyadics.  In  general  it  is  not  possible  to  obtain  the  Green  dyadics  into  closed  form.  Other 
basic  solutions  are  plane  waves.  Sometimes  an  electromagnetic  field  problem  can  be  simplified 
by  decomposing  the  fields  in  two  components.  Each  of  these  components  then  propagate  in  a 
’’simpler”  medium  for  which  the  Green  dyadics  are  known.  Another  way  to  solve  field  problems 
is  the  use  of  duality  transformations  which  allow  us  to  transform  field  solutions  in  one  medium  to 
those  in  another  medium.  Our  aim  was  to  find  the  most  general  media  for  which  decomposition 
of  the  fields  was  possible,  for  which  we  could  find  the  plane  wave  solutions  and  for  which  the 
Green  dyadics  could  be  obtained  in  closed  form.  It  turned  out  that  these  three  problems  are 
intimately  interrelated.  The  common  backbone  behind  these  problems  is  the  possibility  to 
factorize  the  fourth  order  ’’Helmholtz  determinant  operator”. 

This  study  resulted  in  a hierarchical  ordening  - which  we  call  a pedigree  - of  ever  more 
general  bianisotropic  media.  Basically  this  pedigree  consists  of  two  separate  branches.  In  the 
present  contribution  we  will  first  discuss  some  basics  such  as  Green  dyadics,  factorization,  plane 
waves,  decomposition  and  duality.  Then  we  will  focus  on  the  the  two  branches  of  the  pedigree. 

This  paper  gives  only  a small  overview  of  the  subject,  a more  rigorous  historical  overview 
with  more  references  can  be  found  in  [1]. 
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2.  Green  Dyadics  and  Factorization 

The  electric  Green  dyadic  Gge( r)  is  defined  as  the  relation  between  the  electric  current  density 
J(r)  and  the  electric  field  E(r): 

E(r)  = -juj HI Gge(r  - r')  • J(r')dF'.  (2) 

V 

In  a bianisotropic  medium  this  Green  dyadic  satisfies  the  equation 

S£e(V)-Gee(r)  = -I5(r),  (3) 

with  ^ce(V)  the  vector  Helmholtz  operator  given  by  [2] 

£ee(V)  = -(V  x I - ju£)  ■ if1  • (V  x l + juQ  + w2e.  (4) 

To  solve  (3)  it  suffices  to  find  a scalar  Green  function  G(r)  that  is  solution  of 

det^e(V)G(r)  = -«5(r),  (5) 

with  detiJef  (V)  the  Helmholtz  determinant  operator.  The  solution  of  (3)  then  follows  from 

Gee(r)=[^e(V)]TG(r),  (6) 

with  HAee(V)  the  adjoint  operator  of  the  vector  Helmholtz  operator.  The  Helmholtz  determinant 
operator  turns  out  to  be  a fourth  order  operator,  which  makes,  in  general,  a closed  form  solution 
of  (5)  impossible.  For  some  classes  of  media  it  is  possible  to  factorize  this  operator  as  a product 
of  two  second  order  operators,  i.e. 

det  flee(V)  = Ha(V)Hb(V).  (7) 

A medium  for  which  this  is  possible  is  called  factorizable.  Factorizability  does  not  necessarily 
mean  that  we  can  solve  (5)  in  closed  form.  However,  often  a closed  form  solution  is  possible  or 
an  elegant  series  or  integral  representation. 

Sometimes  one  can  factorize  the  second  order  dyadic  Helmholtz  operator: 

&ee(V)=Ha(V)-Hi(V)  (8) 

where  Ha(V)  and  HbCV)  are  first  order  dyadic  operators.  For  these  media  it  is  possible  to  write 
Gee  in  closed  form. 

3.  Plane  Waves 

For  a plane  wave  of  the  form  E(r)  = Eo  exp(— jk  • r)  the  vector  Eo  satisfies  the  equation 

H^e(-jk) -E0  = 0.  (9) 

A solution  different  from  zero  is  only  possible  when  k satisfies  the  dispersion  equation 

detiU-:/k)  = o.  (10) 

If  we  write  k as  ku  with  u a unit  vector  defining  the  phase  velocity  propagation  direction  then 
(10)  is  a fourth  order  polynomial  equation  in  k and  the  dispersion  surface  will  be  a fourth  order 
surface.  When  detH^  is  factorizable  the  dispersion  surface  consists  of  two  second  order  surfaces, 
i.e.  of  two  quadrics.  For  a given  value  of  k the  solution  of  equation  (9)  gives  the  polarization  of 
the  plane  waves. 
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4.  Decomposition 

With  decomposition  we  mean  that  the  fields  can  be  split  in  two  components  as 

E = Ea  + Eft,  H = Ha  + H6.  (11) 

Both  a and  b are  solution  of  Maxwell  equations  for  simpler  media  than  the  original  medium. 

These  simpler  media  are  called  the  equivalent  media.  The  most  well  known  field  decomposition  is 
the  TE  and  TM  decomposition  for  uniaxial  anisotropic  media  [3]  for  which  the  equivalent  media 
are  isotropic  media.  This  can  be  generalised  to  media  where  the  decomposed  fields  satisfy  the 
conditions 

ai  • Ea  + a2  • Ha  =0,  bi  • Ej,  + b2  • H&  = 0,  (12) 

where  ai,  a2,  bi  and  b2  are  arbitrary  vectors. 

Another  way  to  decompose  the  fields  is  the  Bohren  decomposition  [4].  In  this  case  the  fields 
are  decomposed  as 

E = E+  + E_,  H = H+  + H_,  (13) 

with 

E±  = y±H±)  (14) 

with  Y±  some  scalar  constants.  The  most  well  known  medium  that  allows  such  a decomposition 
is  the  isotropic  chiral  medium.  Also  in  this  case  there  are  much  more  general  media  that  allow 
a Bohren  decomposition. 


5.  Duality 

A duality  transformation  transforms  original  fields  E and  H into  dual  fields  E^  and  H<*  as 


E „ 
Hd 


A B 
C D 


(15) 


with  A,  B,  C and  D arbitrary  constants.  The  dual  fields  satisfy  Maxwell  equations  in  a dual 
medium  defined  by  its  medium  dyadics  as 


( & id 

\ Ld  (£d 
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(16) 


Using  the  duality  transformation  it  is  possible  to  transform  field  solutions  from  the  original 
medium  to  the  dual  medium  or  vice  versa.  For  example  the  Green  dyadics  between  both  media 
are  related  through  [5] 


f G-ee,d  G-em,d 
l G-me,d  —mm,d 


A B 
C D 


G*e  G,m 
c c 

—me  1 =Lmm 


0 1 
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(17) 


When  a medium  is  invariant  under  a duality  transformation  we  say  it  is  a self-dual  medium. 
Only  a restricted  class  of  bianisotropic  media  are  self-dual. 


6.  The  First  Branch 

It  turns  out  that  the  class  of  media  that  are  self-dual  or  that  allow  a Bohren  decomposition  or 
that  allow  a factorization  of  the  vector  Helmholtz  operator  is  one  and  the  same.  These  media 
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have  been  studied  in  [6]— [11].  The  isotropic  chiral  medium  is  the  simpliest  representative  of  this 
class.  The  most  general  medium  of  this  class  is  given  by 

e = ea,  p = pa,  £ = (X°L  + m)Vw,  £ = (x«  ~ (18) 

with  a and  k arbitrary  dyadics. 

The  Green  dyadics  for  this  medimn  can  be  written  in  closed  form  [11]  as  follows 

jr)  exp  (fca+  • r)  fexp{-jkD+)\ 

£"(r)  = i+(V)  l J 


jj] exp  (fca_  • r)  /exp  (~jkD_)\ 

2cos e —K  ’ \ Ai\kD-  )' 

(19) 

with 

k = Uy/eji,  = D±  = ,JdetS±^r- ST±  -r,  cos0=^l-x2, 

(20) 

4t(V)  = VV  ± fc(V  • S±)  x 7 + k2Sg1d& tS± 

(21) 

and  where  S±  and  a±  follow  from  a decomposition  of  cos  da  ± k in  a symmetric  and  an  asym- 
metric part  of  the  form 

cos  9a±K  = S_±+a±  x 7.  (22) 

For  more  information  on  the  media  in  this  branch  we  also  refer  to  another  paper  in  these 
proceedings  [12]. 

7.  The  Second  Branch 

The  second  branch  is  much  more  complicated.  This  branch  contains  the  media  that  are  decom- 
posable. It  turns  out  that  these  media  are  also  factorizable.  It  took  us  a long  while  to  find 
the  medium  that  has  all  decomposable  and  factorisable  media  as  special  cases.  This  medium  is 


described  by  rather  complicated  medimn  parameters  [13] 

e = a(z  x I + ajbi  + bjai)  + tj(-Bt  + a2b!  + b2ax),  (23) 

£ = »?(x  X L + a2b2  + b2a2)  + a(B  + aib2  + bia2),  (24) 

£ = r(z  x 7 + aibi  + biax)  - a(-BT  + a2bx  + b2ai),  (25) 

p = — a(x  x 7 + a2b2  + b2a2)  + t(B  + aib2  + bia2),  (26) 


where  a,  r),  r,  ai,  a2,  bi,  b2,  x,  z and  B are  arbitrary  scalars,  vectors  and  dyadics.  It  turns  out 
that  this  medium  is  also  closed  with  respect  to  duality  transformations.  This  means  that  it  is  not 
possible  to  generalize  this  medium  further  using  a duality  transformation.  These  observations 
made  us  conclude  that  (23)-(26)  is  the  most  general  medium  that  allows  decomposition  and 
factorisation. 

Special  cases  of  the  medium  (23)-(26)  were  studied  in  [14]-[34].  An  important  special  case 
are  the  anisotropic  media.  The  most  general  anisotropic  medium  that  allows  factorization  and 
decomposition  is  e = -rpT  + ab  [24],  with  as  special  case  the  uniaxial  anisotropic  medium. 

Another  interesting  class  of  media  are  the  equivalent  media  of  the  medimn  (23)-(26).  These 
are  of  the  form  [13] 

e = —r)BT  + a(z  x 7),  p = tB  - a(x  x 7),  (27) 

£ = aB  + r](x  x 7),  £ = <*RT  + r(z  x 7).  (28) 
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Without  loss  of  generality  we  assume  that  tt/  + a2  = 1.  The  Helmholtz  determinant  operator 
for  these  media  can  be  written  as  a square.  This  allows  a closed  form  Green  dyadic  given  by 

Gee(r)  = {£[(aV  xL-  juB)  ■ (rV  - jux)]  x I 

JUJ 

__L[(^V  + ju;z)(rV  - Jujx)  + a2VV  + jujaV{BTxI)  + 1 x (V  • B)  - (jJ2(B~l)TdetB]} 

JLO 

„ju;c.-R~ 1 -T  „~jkD 

e \ (29) 

ju^i/detRy/D' 

with 

R=  \{a2  + rr]){B  + EF),  c = \[r)B  •x  + rz-H-azxx  + a(detB){B~1)TxI],  (30) 
Z z 

and  

k — cj\J z • H • x + det 7?  — c • R-1  • c,  £)  = ^/r  • R~l  ~ r,  (31) 

where  A*  I is  shorthand  for  xi-Uj  + Uj  xi-us  + u2  x A ■ u2.  The  fact  that  the  Helmholtz 
determinant  operator  is  a square  also  means  that  the  two  second  order  dispersion  surfaces 
coincide.  Each  of  the  equivalent  media  of  a certain  original  medium  have  as  coinciding  dispersion 
surfaces  one  of  the  two  dispersion  surfaces  of  the  original  medium. 
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